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1. Introduction 

The breaking of the strong isospin symmetry exists widely in hadronic physics. For example, 
the charged pion/kaon and neutral pion/kaon form a strong isospin multiplet. Their masses would 
be identical if strong isospin symmetry is preserved exactly. But we observe m 2 n+ — m 2 ^ ~ m 2 K+ — 
m K° ^ 0> which means this symmetry is broken in nature. The mass splitting in the hadron spectrum 
is due to different mass and charge of the fundamental quarks which compose the pion(kaon). 
The mass splitting in the hadron system can be used as a signature to detect the effect of QED 
interactions in QCD phenomena and has been studied by the lattice QCD community in recent 
years [fTL 0, E$J|. It also provides a way to determine the non-degenerate u, d quark masses, which 
can not be made from numerical simulations including QCD interaction alone. The mass of the 
lightest quark, as a fundamental parameter of the Standard Model, is related to the strong CP 
problem and other important questions in high energy physics. 

The mass of the pseudoscalar meson is determined by both QCD and QED dynamics. It is 
difficult to calculate the pseudoscalar mass spectrum theoretically, because of the strong coupling 
in the QCD sector, which makes perturbation theory invalid. Lattice QCD, on the other hand, 
is suitable tool to study hadronic physics non-perturbatively. We employ 2+1 flavor configura- 
tions generated by RBC/UKQCD collaborations [Q, ||, ^] in this work, which include not only two 
degenerate light quarks, but also a heavy strange quark in the sea quark sector. Current lattice 
simulations can determine the pseudoscalar meson mass to roughly the 1% level(statistical error 
only). With the help of the average of ±e trick [f7|], we can determine the mass-squared difference, 
8m 2 = m 2 (QCD + QED) — m 2 (QCD), to even higher accuracy. All of these enable us to build 
the connection between numerical simulation and theoretical description. For example, our result 
can be used to fit the Low Energy Constants (LEC's) in QCD+QED Partially Quenched Chiral 
Perturbation Theory (PQ^PT), which are difficult to obtain from experiment alone. 

Preliminary results have been reported in [||, ^J. In this work, we present a new treatment of 
the violation of chiral symmetry and its further influence on the LEC's and quark mass. We first 
give the PQ^PT description of the pseudoscalar meson mass in Sec. ||[ Explicit chiral symmetry 
breaking in the lattice regularization is discussed in Sec. H The fit of QED LEC's and masses are 
shown in Sec. |j. The summary and plan for future work is given in Sec. ^ l 

2. QCD+QED Partially Quenched Chiral Perturbation Theory 

Chiral Perturbation Theory(^PT) is a low energy effective theory of QCD based on chiral sym- 
metry. QCD+QED Partially Quenched Chiral Perturbation Theory(PQ^PT) gives the contribution 
of the sea quark and valance quark to the pseudoscalar meson mass separately, which gives us more 
freedom to choose the mass of the sea and valance quarks in our lattice simulation. Suppose we 
have a kaon-like meson which is composed of u and s quarks. The mass and the charge of the 
quarks are m\, q\ and m^, qj respectively, = 2e/3, —e/3 for up and down/strange quark. So the 
square of the mass is [[H]] : 

2 2C 2 

m 2 =X u + - I q 2 n + — T (2.1) 
^0 P 

'The mass splitting in the nucleon system will be reported separately in the future. 
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where %n =B Q {m\ +m 3 ), qn = (q\ 
5( 4 ) 



F 2 



(48L5-24l5)to 



-93), and 
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+A(Zl3)<7?3 +4fl(^ y ,Zl3, ^13)913^13 -4fil(Z r ,*13,Zl3)<7?3*13 
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Xl3 



1*1 



(2.2) 



The definition of symbols used in formula d22h can be found in [pi]]. C is the 6{e 2 ) low energy 



constant (LEC). Our qi is eg; of [10]. Our 0{e m) LEC, F,, (/ = 1,2, • •■ ,5), are written in terms of 



linear combinations of Kj of 1 1C]. 



3. Violation of chiral symmetry on the lattice 

Domain wall fermions (DWF) live in five-dimensional space and preserve exact chiral symme- 
try when the size of the extra dimension, L s , goes to infinity. But the size of L s is always finite due 
to the limitation of the computational cost. For example, we set L s =16 and 32 in our work, which 
means we do not have exact chiral symmetry in our simulation. The violation of chiral symmetry is 
quantified by a small additive shift to the input quark mass, the residual mass, denoted as m res [§]. 
Thus the chiral limit is defined as ray + ra res = 0, which leads us to replace all of the quark masses 
in the expression for the pseudoscalar mass in PQ^PT by ra,- + ra res 

The analysis above is based on QCD interactions. If QED is included, it works as a pertur- 
bation to the residual mass which shifts ra mres by ff{e 2 ) from its QCD value. Following usual 
perturbative renormalization arguments, the leading QED effect can be written as: 

Ara res ee ra res (QCD + QED) - ra res (QCD) = C 2 (q\ + q\), (3.1) 

where q\ and q^ are the charges of the quarks composing the pseudoscalar meson. The value of C2 
depends on L s and vanishes exponentially as L s goes to infinity. Higher order terms can be ignored 
at the order we are working, G(e 2 m). 

We next consider the effect of the residual mass to the pseudoscalar meson mass in the 
QCD+QED simulation. The QED contribution to ra res comes into the mass-squared of the pseu- 
doscalar meson when combined with the leading order QCD term, 2BoC2{q 2 + q 2 )- We have two 
different strategies to count this effect. One is to measure the residual mass in the usual way on 
both QCD+QED and pure QCD configurations from the Ward-Takahashi Identity to get Ara res , and 
extract C%. Then we use 2BoC2{q 2 + q^) to account for the QED induced violation of chiral sym- 
metry. The other way is to add a term d jmes (q 2 l + q^) to the expression of ra 2 of the pseudoscalar 
meson, which means we leave 2BqC 2 as new free parameter in PQ^PT formula. We can compare 
whether 8 mres is consistent with 2BqC 2 , which is a way to check whether our understanding of the 
QED violation of chiral symmetry is correct or not. 

We test our understanding about the violation of the chiral symmetry on 16 3 lattice config- 
urations with L S =16 and 32. The result is shown in figure The upper line corresponds to the 
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16 J Ls=16 and 32 result, fit range 0.01-0.02 
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Figure 1: 5m 2 of dd meson with Ls=16 and 32. dmres(...) means 5 mre s(<?i + q\)- It's value at m ps = 
comes from extrapolation. 2B0C2(...) means 2B(,C2{q\ + q\)- (It is shifted horizontally towards left to make 



plot clear.) C2 is determined from the change of the residual mass by setting e ^ and e ■ 
of this term comes mostly from the error in Bq. 



0. The error bar 



L S =16 case; we measure 8m 2 = m 2 (QCD + QED) — m 2 (QCD) of the dd meson with m va i =0.01, 
0.02 and 0.03 and physical charge and extrapolate tom/ + ra res (QCD) = to get 8 mres (q 2 + q\ )■ 
We also determine C2 by measuring the change in m res when e / and e = 0. Then the value of 
2BoC2(q 2 + q\) is calculated and shown on the plot (shifted left a little to make the plot clear). The 
error of 2BoC2(q 2 + ^3) term mostly comes from the error on Bo- It is clear from figure p] that these 
two results are consistent with each other. This consistency is also checked for the L s = 32 case. 
Because the violation of the chiral symmetry of DWF is decreased when L s goes to larger value, it 
is expected that the contribution of the 8 l - mes {q\ + q 2 ) and 2BoC2(q 2 + ^3) should be suppressed at 
larger L s . It can be seen from the plot that both the and C2 terms are decreased by about an 
order of magnitude at the larger value of L s , but their contributions are still consistent within the 
error bar. All of these confirm the validity of using Sfnjes^i ~\~ q 2 *) to account for the contribution 
of the violation of the chiral symmetry. 

4. Numerical Result on Low Energy Constant and Quark Mass 



Formula (2.1) can be used to write the mass-squared splitting as: 



8m = 8m 2 (physical) + 5m 2 (lattice artifact) 



where: 



(4.1) 



W(physical) = 2Cqj 3 /F^ +A(xu)gv3 +^B(Xy,Xi3,Xn)q 2 i3Xi3-^B l (Xr,Xi3,Xi3)q 2 nXi3 
+C[-48L^ 2 3 ^i - \6L r 5 qj 3 Xi3+2A(xis)qi s qi3-2A(x3 S )q3 S qi3]/F^ 
-Yi4Q 2 %i3 + Y 2 4q 2 p x p + F 3 4^ 2 3^i3 - Y^q^Xn + Y 5 \2q 2 n Xi (4.2) 

and 8m 2 (lattice artifact) = 5 mres (q 2 + q 2 ) 
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By fitting the results of 8m 2 , we can obtain all of the QED LEC's, including 5 mres . The 
fit results are shown in figure ||| along with the unitary data points (though all of the partially 
quenched data points within mass cut were used in the fit). From this figure, we can see that 
charged meson receives more finite volume corrections than the neutral meson, which is consistent 
with the theoretical investigation pip. We perform fits with and without the log terms to probe the 
effect of the chiral logarithms. % 2 jdof is adequate in either case, except for 24 3 when we include 
quark masses larger than 0.02 in the fit. One of the reason is that the % 2 jdof for the pseudoscalar 
meson mass fit is around 2 at heavier quark mass, which is a little bit large. The other reason is that 
the quark mass m q =0.02 or 0.03 may be too heavy for SU(3) PQ^PT [Q]. We need to choose the 
fit range carefully. If the heavier sea quark points are omitted, % 2 becomes reasonable. The main 
effect of the logs is to significantly reduce the value of the charged meson splitting in the chiral limit 
(so-called Dashen term). Note, the log terms vanish for the neutral meson splittings at NLO. We 
put the best result of QED LEC's as we have in Table [j] and compare with the phenomenological 



value [10]. Our LEC's come from 24 3 lattices fit , which receive less finite volume correction. The 
fit range is 0.001-0.01 2 , which make SU(3) PQ^PT to be valid. The log terms are also included in 
the fit. 





10 6 C 


10 2 F 2 


io 3 t 3 


io 3 y 4 


10 2 F 5 


1 5 mre s 


X 2 /dof 


24 3 fit 
Ref. ||Tg] 


0.27(19) 
7.3 


1.59(10) 
0.38 


-10.6(7) 
1.58 


9.8(16) 
2.83 


2.00(68) 
-0.953 


5.08(9) 
N/A 


2.11(73) 
N/A 



Table 1: The preliminary result of the QED LEC's from fits of 8m 2 and only statistical error is quoted. The 
QCD LEC's were taken from the SU(3) PQ^PT fit by RBC/UKQCD collaborations. All of the QCD and 
QED LEC's are defined at scale jX = lGeV. The last row gives phenomenological estimate in ][Io|]. 

Next, we determine the physical u,d, and s quark masses from the QCD LEC's given by the 
RBC/UKQCD collaborations and the QED LEC's from our fit. The fluctuation in QCD configu- 
rations is considered by including the jack knife blocks of QCD LEC's. Since there are three un- 
known quark masses, we need the experimental value of the mass from three pseudoscalar mesons 
(7T + , K + and ^°) to fix the quark masses [13]. The mass of the neutral pion is not used, because 
we don't include the disconnected diagram in our simulation due to the difficulty of computation. 
The renormalized quark masses are defined as m^ s = Z m (m q + m res ), where Z m is the matching 
factor from lattice to MS(NDR) scheme at a scale /I = 2 GeV obtained via the RI/MOM scheme 
with the non-perturbative technique 3 . We use the value Z m = 1.656(48)(150) from [@]. Table § 
gives the result of the quark masses. The analysis on the 24 3 lattice has less finite volume effect and 
reasonable % 2 for the LEC's fit for the fit range 0.001-0.01. So we quote the quark mass determined 
from these parameters as our preliminary result, m™ s = 2.79(37)MeV, mf s = 4.84(52)MeV and 
mp = 95.9(9.6)MeV. 



5. Summary and Plan of Future Work 

Our simulations are based on 2+1 flavor QCD configurations generated by RBC/UKQCD 

2 We add m va ] = 0.001 data to our analysis. So our fit result are changed form those reported in ^ as well the result 
quoted in fel], where only 0.005 and heavier quark masses are used. 



Here our Z m is calculated on pure QCD configurations. 



5 



Isospin breaking in 2+1 flavor QCD+QED 



Ran Zhou 



16 lat. fit range:0.01-0.03 



16 lat. fit range :0. 01-0.03 





24 lat. fit range:0.001-0.03 





Figure 2: 16 3 linear fit(upper left panel), 16 3 chiral log fit(upper right panel), 24 3 linear fit(lower left panel) 
and 24 3 chiral log fit(lower right panel) lattice data points and fit results for the meson mass splitting. Lines 
correspond to fits with and without chiral logs. The data points in the plot correspond to ud,uii and dd 
mesons, respectively, from top to bottom. 



lattice 


fit range 


m u 


m d 


m s 


m„ 
m d 


m s 
mud 


24 3 


0.001-0.01 


2.79(37) 


4.84(52) 


95.9(9.6) 


0.57(1) 


25.1(5) 



Table 2: Preliminary results of the u, d and s quark masses determined from pure QCD interaction. The 
value is in the unit of MeV and MS(NDR) scheme at renormalization scale fi = 2 GeV. The QCD LEC's 
were taken from the SU(3) fit by RBC/UKQCD collaborations and QED LEC's are from Table 0. The error 
quoted here is only statistical error. 



collaborations and quenched non-compact QED configurations generated by us. After treating the 
violation of chiral symmetry carefully, we fit the pseudoscalar meson mass-squared differences 
using QCD+QED PQ^PT to extract the QED LEC's. The LEC's are affected by the volume of 
the lattice, fit range and the formula used in the fit. We have investigated these systematics, but 
have not yet quoted errors for them in our preliminary results. The up, down, and strange quark 
masses are also determined using the LEC with the experimental mass values for pseudoscalar 
mesons (pi+, K+, K0). Our preliminary results for MS(NDR) scheme at the renormalization scale 
jX = 2 GeV is obtained using LEC's determined on 24 3 lattice and fit range as 0.001-0.01. They are 
m]p = 2.79(37)MeV, mf* = 4.84(52)MeV and mf* = 95.9(9.6)MeV. Errors are statistical only. 
In this work, quenched QED configurations were used to account for the EM interaction. The 
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systematic error caused by this approximation can be removed by the reweighing method [|12|]. In 
addition, the analysis on different lattice spacing will allow to extrapolate to the continuum limit. 
This will be done on the companion 32 3 lattice ensemble generated by the RBC and UKQCD 
collaborations (see the talks by Kelly and Mawhinney in these proceedings). 
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